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Abstract — An observer-based Hamiltonian identification al- 
gorithm for quantum systems has been proposed in [2]. The 
later paper provided a method to estimate the dipole moment 
matrix of a quantum system requiring the measurement of 
the populations on all states, which could be experimentally 
difficult to achieve. We propose here an extension to a 3-level 
quantum system, having access to the population of the ground 
state only. By a more adapted choice of the control field, we will 
show that a continuous measurement of this observable, alone, 
is enough to identify the field coupling parameters (dipole 
moment). 

Keywords: Nonlinear systems, quantum systems, pa- 
rameter estimation, nonlinear observers, averaging. 

I. Introduction 

For applications ranging from quantum computers to the 
synthesis of new molecules, an accurate estimation of the 
parameters involved in the dynamics of the quantum system 
is fundamental. Various methods have been engineered 
over the years Uke the maximum-likelihood methods [10], 
[8], [11], the maximum-entropy methods [4] and minimum 
Kullback entropy methods [13]. The optimal identification 
techniques via least-square criterias [6], [5], [12] and the 
map inversion techniques [14] are some other techniques 
explored in this area. In [9], a state-observer is presented 
for the state identification combined with a gradient method 
on the dipole moment. This result was then improved in 
[2] succeeding in simultaneously estimating the state of 
the system and it's dipole moment using observers. In [12] 
a rigorous proof of the well-posedness of the problem is 
proposed. All these results required the knowledge of the 
populations on all energy levels. Experimentally, that is 
extremely difficult to achieve. Since in quantum mechanics, 
measuring an observable influences the system, the less 
information we need, the less we disturb the system, and 
the more likely our estimation is accurate. Our goal was to 
improve the result given in [2] in order to estimate the dipole 
moments of a quantum system measuring continuously the 
population on the first state only. We focus here on 3-level 
systems with a single population measure. 
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In section [III we explain the 3-level system and set the 
estimation problem attached to ([T]i. In section |lll] we present 
a 2-step estimation procedure based on two nonlinear 
asymptotic observers ^ and (|4]i. Sections |IV] and |V] are 
devoted to local convergence proofs. 

II. The 3-level system 
A. Model and problem setting 

Denote by \k), k= 1,2,3 the 3 states of energies such 
that |£'2 ^ £^1 1 7^ 1^3 — ^2 1 ■ Throughout the paper we use the 
following notations for k,l = 1,2,3: a" = |/) - \k) {l\, 
ai' = |/} (^1 + \k) {l\, (jf = |/} (/| - \k) {k\ and P, = |^) {k\ 
(projector on \k)). Assume that the dynamics is described 
by the following Schrodinger equation: 

1 1*?) = ^ (Ho + A(f )//,) 1^) , {^'\P, \^') , 

where |^) is the wave-function, A(f ) e R the electromag- 
netic field, Hq — Y}k=i^kPk the free Hamiltonian, H\ = 
fJ-n^^x^ + l^23<^x^ the Hamiltonian matrix describing the 
coupling with the electromagnetic field (dipole moment) 
and y the measurement output. Assuming the energies Ek 
known, the goal consists in estimating the real coupling 
parameters jj.12 and from the output y. We assume the 
electromagnetic field resonant with transitions 1—2 and 
2-3: 

fE2—E\ \ - //is— £2 \ 
A(f) = Mi2Ai2sin f — - — 1\ +M23A23sin f — - — 1\ 

with small amplitude magnitudes A 12 and A23 and nor- 
malized slow modulations |mi2|,|m23| G [0,1]. We have 
|Ai2MpM^'3M23| ^ \E2- E\\,\E^ - E2\. In the interaction 
frame |<I>) = eT"' and after neglecting highly oscillating 
terms (rotating wave approximation) we get the following 
model 

^ |4>) = (Ml2ni2(7^2 + M23n23CT23) \^) 

where ^12 = and ^23 = are Rabi ampUtudes 

when (mi2,M23) = (1,0) and (mi2,M23) = (0,1). 

In the sequel we will use the density operator p = |<I>) (<I>| 
instead of the wave function |<I>}. The estimation of the real 
parameters /X12 and /X23 is then equivalent to estimation of 
the two other real parameters £l\2 and Q.23 appearing in the 
dynamics of the projector p 

^P =Mi2i2l2 [C^'^P] +M23ii23 [(^^^P] (1) 



via the output y = Tr(Pip) and using ui2 and M23 as 
excitation real inputs. Remember that cr'^ and a^^ are anti- 
symmetric and real matrices: if the entries of p are initially 
real, they remain real; if p is initially a projector and thus 
describes a pure quantum state, its remains a projector. 
Since we are in the 3-level case, p can be seen as a point 
on the two dimensional manifold MP^, the projective space. 

B. Identifiability 

It is proved in [12] that it is possible to identify £212 
and £^23 by measuring all the populations, i.e., via the 
measurement outputs (Tr(fip) ,Tr(P2P)) (Tr(P3p) = 1 — 
Tr (f 1 p ) - Tr (f 2P ))■ With just y = Tr (Pi p ) , we provide here 
below arguments showing identifiability via adapted choices 
for inputs M12 and M23. With ui2 = l and M23 = we recover 
essentially a 2-level system with states |1) and |2) and we 
can identify £2i2 from y following [2]. This corresponds to 
the first step that is treated in theorem [T] 

Assume now that £2i2 is known. Set M12 — 1 and M23 — 
rj cos 9 with -^9 — £ii2 and 77 a small positive parameter 
With ^ = e^^^ pe^*^ , the output map becomes 

_ Tr((Pi +/>2)^)+cos(2e)Tr((T^2^) + sin(2e)Tr (gj^^) 

and obeys to 

= ^co^9 £223[cT^^^] + I?cos0sin0 ^223[ff^^^] 
Since 77 <C 1 we can average its dynamics: 

'-^A^a-\^\. (2) 

The average values of 3'(f)(l + 2cos(20)) and y(t)(\ — 
2 cos (20)) are Tr(Pi<^) and Tr(P2^), respectively. Thus in 
average all the populations are measured and according 
to [12], £^23 is identifiable. This second step is treated in 
theorem |2] 

III. Estimation algorithms and simulations 

As explained here above, we proceed in two step. In a 
first step we set in ([U, M12 — 1 and M23 — and estimate 
from the output y(t) the parameter f2i2 via the following 
nonlinear dynamical system (an invariant nonlinear observer 
inspired by [2], [1], [3]): 

^p=£ll2[fTl2,p] + ... 

...eri2()'(f) -Tr(/'ip))((T^2^ +p(T^2 _2Tr {o]}p) p) 

|:n,2 = e2yi2Tr(cTi2[c,i2,p]) _Tr(Pip)) (3) 

with ri2,7i2 positive parameters of order 1 and £ a 
small positive parameter. Local convergence is proved in 
theorem[T] The dynamics ([3]) respect two important features: 
if the entries p are initially real, they remain real for f > 0; if 
p is initially a projector and thus describes a pure quantum 
state, it remains a projector for f > 0. 

Assuming £2i2 obtained via this first step, we take, as 
explained in previous section, M12 = 1 with M23 = cos 9 



(^0 = £^12 and 77 a small positive parameter) to estimate 
£223 via a second nonlinear dynamical system 

^p = £ll2[fT'lpl + 77COse n23[(72^pl + ... 
at 

...£T]r23(3'-y)(l-2cos(20))(zfp + pZf -2Tr(z23p)p) 
^^23 = e2r7r23(.v-j)(l-2cos(2e))Tr(zf [Z23,p]) (4) 

where Y?'^ = V (t)cf-'^U\t), Ef = U[t)aV'U\t) with 
U{t) = exp(0(f)(j'^) and where r23 and 723 are positive 
parameters of order 1 and e is a small positive parameter. 
Local convergence is addressed in theorem ID 
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Fig. I. Estimation of i2i2 in a first step (t e [0,50]) via (5); estimation 
of Q.23 in a second step (t > 50) via (no modeling and measure errors, 
il[2 and stand for Q,i2 and ^23). 
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Fig. 2. Similar simulations to those of figure[T]but with 20% of Gaussian 
additive noise on the output y(t) and 10% of Gaussian additive noise on 
the inputs »i2 and M23. 

Let us look at some simulations (figures [T] and |2]) with the 
following numerical values: p(0) = p(0) = Pi, £2i2 = 1-0, 
£223 = 0.8, £2i2(0) = ^, £223(0) = L5£223, e = ?? = |, 



12 — 



= 4, 7i2 = 723 = 1- Convergence of Q.12 is 



with e£ii2 = ^ii2 ~^i2- Since 4-^ = 0, for e small enough 



effective after t > 20 that corresponds to i Rabi periods we can consider the average system: 



associated to transition 1—2 since w 18. Convergence 
of £223 is achieved for t G [50,200]. The interval length 
corresponds to Jr- Rabi periods associated to transition 
2 — 3 for the average dynamics dU since ^^^^ ~ 140. 
These convergence times are in good agreement with the 
convergence times that can be obtained from the linearized 
system Q appearing during the proof of theorem [T] When 
additive noises are introduced on the inputs and output, 
the performance are not dramatically changed and the 
convergence times are almost the same. 

IV. Estimation of 

Theorem 1: Take system ([T]) with inputs M12 = 1, M23 = 
and consider the estimation of p and Q.12 via (O. Take 
ri2,7i2 > and assume that p(0) and p(0) are real projec- 
tors with Tr((Pi +P2)p(0)) e]0,l[. Then, for e > small 
enough, exists a > such that, if 1 - Tr(j5(0)p(0)) < 
(7 and |ni2(0) - £2i2| < C7, then limrH^+«,p(f) - p(f) = 
and lim,H^+oo£2i2(0 = Qii- Moreover the convergence is 
exponential. 

Proof: Since p and p remain real projectors for f > 0, 
they can be seen as points on the two dimensional manifold 
RP^, a projective space. In particular ([T) is a dynamical 
system with only 2 degrees of freedom whereas the space 
of 3 X 3 symmetric real matrices where the dynamics is 
expressed is of dimension 9. We have try to use less 
scalar variables but averaging computations performed here 
below are then much more complicated. In fact calculations 
based on 3 x 3 symmetric real matrices and thus with more 
variables than necessary simplify notably the analysis. 
Set ^9 = £2i2 and consider the unitary and real transfor- 



mation (P3 = |3) (3|) 
f/(f) =exp(0(7'2) 



P3 + cos e (Pi + P2) + sin 0(7 



12 



and the attached change of frame ^ = pU, £, = pU. 
Since 



U'^PiU 



P1+P2 cos(20) p sin(20) 
— :^ ^ ^ — <^x 



and 



system (|3]l reads 
d 



dt 



P1+P2 cos(2e) 12 

' CT- 



sin(2e) 12 

z CT, 



2 2 
cos(20)d2 + sin(2e)c7^'2^| 

-h| (^cos(20)cTjf2 + sin(20)o;l2 
-2Tr (^(^cos(2e)fT^2_^sin(2e)fT^'2j |j | 

^-f2l2 = eyi2Tr(^(^cos(2e)(Tj^ + sin(20)o-i 



d - 



.Tr 



P1+P2 



cos(26) 19 
^ — -(J, - 



sin(20) 



d s 
dt 



^=eni2[cT'^^] 



e^Tr I)) + Icri^ - 2Tr (cri^l) | 



.e^i^Tr (ct12(^ - |)) (cr^^l + Icrf - 2Tr [o}^) | 



|n,2 = afTr(cTi^[cT'Ml)Tr(cT^^(§-|; 

.712. 



.a^Tr(cTi2[cT'Ml)Tr(cTf(^-|: 



But ^ = ^ and £2i2 = is a steady state of this average 
system. Assume we have proved that this equilibrium is 
exponentially stable. Then the averaging theorem (see, e.g. 
[7, theorem 4.1.1, page 168]) ensures that the above time- 
periodic system admits a unique periodic orbit exponentially 
stable near (^,0). Since (<^,0) is also an equilibrium of 
this time-periodic system, this exponentially stable orbit 
coincides with this equilibrium and the theorem is proved. 

Let us prove now that {£,,0) is a hyperbolically 
stable equilibrium of the average system. We have 

Tr(o;i2[cjiM]) = Tr ([a^\a'^]l) and [a^^a'^] = 
-lal^. Thus Tr(o;i2[(7i2^|]^ = 2Tr (q^'^l^ . Similai-ly 
Tr('ff^'2[(yi2 ^ -2Tr(^(7,'2|y Thus the average system 
reads 



dt 



r. 



Alt 



2Tr 



Tr((Ti2|)Tr(cTf(^-| 



(5) 



Assumption Tr((fi +f2)p(0)) > implies that 
Tr((Pi +P2)<^) > 0. We can choose the initial value of 
such that^Tr(A^) = Tr((Pi +^2)^) > and Tr(P2^) = 0. 
^ and I belong to RP^ and around ^ the variables 
X = Tr^cT^^^I^ and £ = Tr^Pi|^ form local coordinates 

for ^: when £, = x = and £ = a with a G]0, 1[. Some 
standard computations yield to the following linearized 
dynamics: 

d ^ ^ ^ eaTn . 

— X — —2eaLl]2 X 

dt 2 

d ea{l — fl)ri2_ 

~r-? — 7^ z 

dt 2 

^n,2 = (6) 

with i = i and z — z — a. This linearized system is expo- 
nentially stable. 



Remark 1: The stability of the above average system Q 
is more than local. It admits the following Lyapunov func- 
tion: 

^ (iinf + Tr ((7^2 (I - ^ )) V Tr (o;'2 (I „ ^ )^ 2 

Even if theorem [T] is a local stability result, the proposed 
estimator (O should have a large attraction region. This is 
corroborated by simulations of figure [T] 

V. Estimation of 

Theorem 2: Take system ([T]) with inputs u\2 = 1, M23 = 
77COS0 where r\ is constant and ^9 =^12. Consider the 
estimation of p and H^i via (|4|. Take r23 > and 723 > 0. 
Assume p is a real projector with Tr((fi +P2)p(0)) > 0. 
Then for e, 77 positive and small enough, exists (7 > such 
that, if p(0) is a real projector such that 1 — Tr(p(0)p(0)) < 
a and |£l23(0) -il23| < C then lim,H^+«, p (f) - p(f) =0 and 
lim(^+oo £223(0 =^2i- 

Proof: The unitary and real transformation U — e^" 
reads P3 + cos (Pi+P2) + sin0 a Consider the attached 
change of frame ^ = U^'pU, | = f/^pf/. Since U^a^^U = 
COS0 (7^^' — sin0 a^^, t, obeys to 

=Tjcos20 £223 [C7^^ ^]-Tj cos sine £223[CT'^^] 
and (in becomes 

4-1 = 7? cos^e fiasi'T^^,!]-?] cose sine £i23[cf'M] + -- 
at 

...e7?r23^Tr((/i2+cos(2e)cTi2 + sin(2e)cTl2) (^-|))... 

...(l-2cos(2e))(cTf I + |cTf-2Tr((7f 1)1) 
^£223 = e'i7to^Tr((/i2 + cos(2e)cTf + sin(2e)cTi2) -|)) ... 
...(l-2cos(2e))Tr(cTf [CT^M]) 

In average, 

Tr ((/12 +cos(20)(7^2 ^ sin(20)(7^'2) - |)) (1 -2cos(20)) 

is equal to Tr^P2(^— 1)^ After neglecting the highly 
oscillating terms, we obtain: 

^^='7^£^23[a-,^] 

^|=,in23[a-|] + ... 

••.er)r23^Tr [p^M - I)) (erf | +|cTf - 2Tr (af |) |) 
|£223 = e'i)r23^Tr (/'2(§ - I)) Tr (erf [a^M]) 



of the above average system converges locally exponentially 
towards (^,£223). Since (^,£223) = (^,£223) is also solution 
of the original system (01, this implies that, for T] small 
enough, ((^,£223) converges locally towards (^,£223). This 
convergence is exponential. H 

VI. Conclusion 

For a 3-level system ([T]) with only a single population 
measurement we have proposed an algorithm in two steps 
for the estimation of £2 12 and £223. Simulations show 
the robustness to additive noise of this algorithm relying 
on nonlinear asymptotic observers preserving the usual 
symmetries (change of frames). Theorems [T] and |2] ensure 
the local and exponential convergence. We can imagine 
switching periodically between estimation of £2 12 via (O 
and estimation of £223 via (|4|i in order to produce estimations 
of £2i2 and £223 in real-time. 

References 

[1] S. Bonnabel, Ph. Martin, and P. Rouchon. Invariant asymptotic 

observers. IEEE Trans. Automatic Control, 53:2514-2526, 2008. 
[2] S. Bonnabel, M. Mirrahimi, and P. Rouchon. Observer-based Hamil- 

tonian identification for quantum systems, to appear in Automatica. 
[3] S. Bonnabel, M. Mirrahimi, and P. Rouchon. Qubit Hamiltonian 

identification: a symmetry-preserving observer-based approach. In 

17th IFAC World Congress, 2008. 
[4] V. Buzek. Lecture Notes in Physics, volume 649, chapter Quantum 

tomography from incomplete data via MaxEnt principle, pages 1 80- 

234. Springer- Verlag Berlin, 2004. 
[5] J. M. Geremia and H. Rabitz. Optimal identification of hamiltonian 

information by closed-loop laser control of quantum systems. Phys. 

Rev. Lett, 89:263902, 2002. 
[6] J. M. Geremia and H. Rabitz. Optimal hamiltonian identification: 

The synthesis of quantum optimal control and quantum inversion. / 

Chem. Phys, 118(12):5369-5382, 2003. 
[7] J. Guckenheimer and P. Holmes. Nonlinear Oscillations, Dynamical 

Systems and Bifurcations of Vector Fields. Springer, New York, 1983. 
[8] R.L. Kosut, I.Walmsley, Y. Eldar, and H. Rabitz. Quantum state de- 
tector design: Optimal worst-case a posteriori performance. Technical 

report, arXiv: quant-ph/0403150 2004. 
[9] R.L. Kosut and H. Rabitz. Identification of quantum systems. In 

15th IFAC World Congress, 2002. 
[10] R.L. Kosut, H. Rabitz, and I. Walmsley. Maximum likelihood 

identification of quantum systems for control design. In I3th IFAC 

Symposium on System Identification, Rotterdam, Netherlands, 2003. 
[11] R.L. Kosut, 1. Walmsley, and H. Rabitz. Optimal experiment design 

for quantum state and process tomography and hamiltonian parameter 

estimation. Technical report, arXiv:quant-ph/041 1093 2004. 
[12] C. Le Bris, M. Mirrahimi, H. Rabitz, and G. Turinici. Hamiltonian 

identification for quantum systems: well-posedness and numeiical 

approaches. ESAIM: COCV, 13:378-395, 2007. 
[13] S. Olivares and M.G. Paris. Quantum estimation via the minimum 

kullback entropy principle. Phys. Rev. A, 76:042120, 2007. 
[14] N. Shenvi, J.M. Geremia, and H. Rabitz. Nonlinear kinetic parameter 

identification throughmap inversion. / Phys. Chem. A, 106:12315- 

12323, 2002. 



In the time scale rjt instead of t and up to a circular 
permutation (2,3,1) to (1,2,3), we recover (O of theo- 
rem [U We can always choose the initial value of such 
Tr(Pi^) = Tr(P2^) = "^''"^'^^'^^^ > 0. Thus assumptions of 
theorem [T] are satisfied, in particular Tr((P2 +^"3)^ (0)) £ 
]0, 1[. Consequently, for e small enough, (^,£223) solution 



